Centrality dependence of strangeness and (anti)hyperon production at BNL RHIC 
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We evaluate strangeness produced in Au-Au interactions at ^/snn = 200 GeV, as function of 
reaction participant number A, and obtain the relative strange quark content at hadronization. 
Strange baryon and antibaryon rapidity density yields are studied, relative to, and as function of, 
participant number, and produced hadron yields. 

PACS numbers: 25.75.Nq, 24.10.Pa,12.38.Mh 



The experiments at the Relativistic Heavy Ion Col- 
lider (RHIC) study the properties of quark-gluon plasma 
(QGP) matter, a state that was present in the early Uni- 
verse. We explore in this paper the production of baryons 
and antibaryons containing one or more strange quarks. 
Strange antibaryons offer a unique signature of the pres- 
ence of QGP [l| . Their yields probe how QGP turns into 
conventional matter (hadronization). 

At origin of this strange (anti)baryon signature is the 
high abundance of strangeness expected in QGP, com- 
bined with the high matter density at hadronization 
of this entropy rich phase 0. A relatively small con- 
ventional production background enhances the impor- 
tance of (multi)strange antibaryon QGP signature 0. 
We use results of an analysis of the general pattern 
of particle production at central rapidity as function 
of A 4] for RHIC y/s^ = 200 GeV Au-Au reac- 
tions. These result were obtained within the statisti- 
cal hadronization model (SHM) description of the 
PHENIX experimental central rapidity yields dNi/dy, 
Ni — tt^~ , , p and p along with STAR experiment rel- 
ative yields K*(892)/K" and </>/K" 0. 

The hadronization conditions for strangeness are pre- 
sented for the first time here and will be used to pre- 
dict the centrality dependent yields of strange and multi- 
strange baryons and antibaryons. These results are ob- 
tained without any adjustment or revision of Ref. 0. 
The statistical model parameters allowing us to differ- 
entiate the degree of chemical equilibration is the phase 
space occupancy 73. 7^ has for each flavor T of quarks 
a similar role as the "blackness" of a gray-body photon- 
emitter: for 7 < 1 there are fewer particles than ex- 
pected from a black-body source, and for 7 > 1 more. 
The SHM analysis of the final state hadron yields pre- 
sumes thermal equilibrium, and we consider all chemical 
(non-)equilibrium conditions studied in literature: 

1) Chemical non-equilibrium of valence light and 
strange quark pair abundance; this we expect to arise 
in a rapid transformation of the deconfined quark-matter 
into free-streaming hadrons. In such a sudden hadroniza- 
tion, there is no time to re-equilibrate the final state 
yields determined by fragmentation and recombination 



of available QGP partons. Because the (chemically equi- 
librated) phase space in QGP, in general, has a different 
density compared to the hadronic gas (HG), QGP chem- 
ical equilibrium described by quark-matter space occu- 
pancies 7^ = 1, 7^ = 1 implies hadron chemical non- 
equilibrium, with the parameters 7^ ^ 1, 7J 1 7^ 1. 

2) Chemical semi-equilibrium, which may arise for suf- 
ficiently slow hadronization, allowing to re-equilibrate 
the light quark abundances after hadron formation, 
7J 1 — * 1. Since strangeness re-equilibration is a consid- 
erably slower reaction process |9J, it did not have suffi- 
cient time to develop. Alternatively, this chemical semi- 
equilibrium can also arise in reactions not involving a 
phase change, with particle production in hadron sector 
evolving to chemical equilibrium for non-strange hadrons, 
but strange hadrons not having enough time to reach the 
strangeness yield 'absolute' chemical equilibrium. 

3) Full chemical equilibrium, 7^ = 7J 1 = 1, requires a 
complete re-equilibration, or, if there is no phase change, 
equilibration of all types of hadronic particles. 

We assess the viability of QGP strange (anti) baryon 
signatures of dcconfincmcnt by comparing these three 
scenarios as function of A. For very peripheral colli- 
sions, also studied in Ref. 0, there was a significant de- 
viation from chemical equilibrium with the strangeness 
yield being significantly reduced, as measured both, by 
quark occupancy 7^ /7J 1 < 1 , and the (nearly conserved) 
strangeness to entropy ratio s/S ~ 0.019. Moreover, in 
this limit, the data fit employing chemical equilibrium 
Js = Jq = 1 did not yield a confidence level P[%] com- 
parable to the chemical non-equilibrium approach. The 
present analysis includes all results for A > 20 which 
yield P[%] > 60%. With increasing A, for the chemi- 
cal non-equilibrium statistical hadronization model, the 
strangeness occupancy 7^ was rising rapidly, reaching 
for most central collisions the value 7^ = 2.4 . The rel- 
ative yield of strangeness per entropy was found for the 
most central collisions, A ~ 350, at s/S ~ 0.028 ± 0.002, 
which corresponds to the expectations one has for nearly 
chemically equilibrated QGP. 

We show, in the top panel of Fig.Q] the average 
strangeness density of the hadron source at central ra- 
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FIG. 1: Top panel: Strangeness density ds/dV [fm -3 ]; 
middle panel: geometrically reduced hadronization volume 
dV/dy/A 2 ^ [fm -3 ]; and bottom panel: the ratio of light 
quark to strange quark density q/s as function of A. Solid 
lines: chemical non-equilibrium; chain lines: chemical semi- 
equilibrium; and dashed lines: chemical equilibrium. 

pidity, as function of A, at the time of hadronization: 
(ds I dy) I (dV / dy) — ds/dV — ds/dV. For very periph- 
eral collision systems, the central rapidity density is at 
ds/dV = 0.03/fm 3 . As A increases, there is a signif- 
icant increase in the strangeness density at hadroniza- 
tion in the chemical non-equilibrium model, with the 
full non-equilibrium (solid line) approaching a five-fold 
increase. ds/dV doubles in the semi-equilibrium case 
(chain line) and increases by 50% in chemical equilibrium 
model (dashed line). 

This result combines two multiplicative factors: 

1) The strangeness yield per unit of rapidity, ds/dy, 
which, in principle, can be directly estimated from the 
experimental yields of particles which were fitted, and is 
thus model independent. 

2) The statistical volume per rapidity dV/dy, which is a 
model dependent result of the fit to the global particle 
yield pattern. 

We show its behavior in the middle panel of Fig.^ re- 
duced by the initial state geometric scale factor A 2 ! 3 . 
For small A, the chemical non-equilibrium result shows a 
40% reduction of the residual volume with increasing A, 
as would be expected for an increasing stopping power, 
e.g., due to the progressive onset of deconfincment. 



All models show, for large A, an increase of 
(dV/dy)/A 2 / 3 with A, which can be interpreted as being 
due to the grater dynamic expansion prior to hadroniza- 
tion of the larger, initially more compressed system. For 
the chemical non-equilibrium model, (dV / dy) / A 2 / 3 in- 
creases by 180% compared to minimum near A = 25. The 
yet larger increases, by a factor 3.5-4, seen for the the 
chemical semi-equilibrium and equilibrium models, com- 
pensate for the chemical non-equilibrium factor ~ 1.6. 
Some may consider such an increase excessive, in which 
case the result is an indirect evidence for > 1. 

In order to gauge the impact of the strangeness 
hadronization density increase, we have to compare it 
to light quark density, as this effectively eliminated the 
model dependence. To obtain such a q/s ratio, we esti- 
mate the light quarks number, including any originating 
in hadronization fragmentation, from the final state en- 
tropy content. For quarks with thermal mass m q ~ 2T, 
for the classical massive gas the entropy per particle is 
S/N ~ 5. Thus: 
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We reduced the total entropy content by the two (k s = 2) 
fractions s and s, and divided by k(S/N) = 20 to obtain 
the approximate average yield of the k = 4 components 
q = u, u, d and d. 

The chemical non-equilibrium q/s (solid line in bottom 
panel of Fig.^l and semi-equilibrium (chain line) ratios 
q/s are identical within the error margin of this estimate. 
Strangeness is, at small A, suppressed by more than fac- 
tor 2. For A — > 350, a 20% excess of light quark num- 
ber compared to strange quark number remains. This 
can be understood, e.g., as being due to the thermal 
strange quark mass being slightly larger than the light 
quark thermal mass, and both light and strange quarks 
approaching chemical equilibrium yields in the QGP. 
The hadron-phase chemical equilibrium result (dashed 
line) shows no variation with A. The increased yield of 
strangeness is consistent with the rise of the reduced reac- 
tion volume, which is implying an increase in the lifespan 
of the dense phase. 

In the following, we discuss results which include decay 
products of weak decays: we count all pions from Ks de- 
cay, we include 10% Kl decay products, and all hyperon 
decay, i.e., both pions and baryons are accepted. In par- 
ticular, we included in our dA/dy and dA/dy, all weak 

decays of dSr/dy, d"E. + /dy and dfl~/dy, d£l + /dy. These 
idealized assumptions will cause minor calculable correc- 
tion for more realistic experimental acceptances. The 
interested reader can use the SHM implementation Q to 
obtain results for other weak decay conditions. 

We present, in Fig.|3 the predicted production rate 
of dA/dy, dA/dy (top panel), dE~/dy, dE + /dy (bottom 
panel) as function of dhr/dy and dn~/dy at central ra- 
pidity. The dependence on h~~ is indicated in top, and 
on 7r~ in bottom of the figure. We note that the equilib- 
rium model (dashed lines) predicts fewer (anti)hyperons. 
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We omitted the chemical semi-equilibrium lines which 
are in-between the two other models. The available 
experimental results shown are for ^/snn = 130 GeV 
Au-Au rather than here considered 200 GeV reactions, 
dA/dy, dA/dy 10], and dS~/dy, dZ + /dy jUj. Consider- 
ing that at y'sNN = 200 GeV the yields should be higher, 
the chemical nonequilibrium model is more likely to be 
consistent with the 200 GeV data. 
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FIG. 2: (color online) The predicted hyperon (blue) and, in 
same panel directly below, antihyperon (red) yields. The ra- 
pidity yields are shown as function of both dn~/dy (bottom 
scale) and dh~/dy (top scale) for t/snn = 200 GeV Au-Au re- 
actions. Solid lines for chemical non-equilibrium, and dashed 
lines for the chemical equilibrium. The experimental results 
(square symbols for hyperons and triangles for antihyperons) 
are from ^/snn = 130 GeV Au-Au reactions, presented here 
for reader orientation. 

We now turn to discuss the enhancement effect in pro- 
duction of hyperons and antihyperons. A way to study 
this enhancement is to consider the yield of the particle 
per hadron and/or per pion as function of hadron/pion 
yield. We present, in Fig-OJ the full set of hyperon the- 
oretical predictions normalizing these with the ir~ yield. 
We see that the normalized chemical equilibrium yield 
is fiat, while the chemical non-equilibrium yields predict 
enhancement with A. This effect can be seen, in Fig.|2 
as a slight up-bend in the particle yield at high A. 

Such an up-bend is visible in the experimental A and 
A results for ^/snn = 130 GeV Au-Au reactions |lflj . 
however the data was not analyzed to reveal this effect. 
Considering result shown in Fig-EI we suggest that a fit to 



the experimental centrality dependence should be made 
allowing for a non-linear shape, e.g.: 

dN/dy = ai h- + 0.5a 2 {h-) 2 , (2) 

where predicted a 2 is, up to a factor h~/n~ = 1.2, the 
slope seen in Fig.[31 solid line. 
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FIG. 3: (color online) The yields of hyperons (blue) and 
antihyperons (red), A, A (top panel), E~ H (middle panel), 
and Q~, fi bottom normalized with tv~ yield, as function of 
dir~/dy yield. The lines are as in Fig. [5] 

Should a.2 — > 0, this would constitute reliable evidence 
for validity of the chemical equilibrium model: such a re- 
sult relates closely to the flat q/s relative yield, in Fig.^ 
bottom panel, ai — > means that there is no central- 
ity dependence in relative strange and light quark abun- 
dances, and densities, at hadronization. This is expected 
for chemical equilibrium, but would require an extraordi- 
nary coincidence to occur at several centralities for chem- 
ical non-equilibrium. Therefore, this is a test for chemical 
equilibrium, more appropriate than the more common 
study of 7^ — > 1 and 7J 1 — > 1. In the latter case, the 
fitted value of yy~„ can be different from unity, due to 
incompleteness of the SHM set of hadron states. 

There is another way to study enhanced production 
of strange hyperons 12j: we consider the hyperon yield 
normalized by A, as function of A. To obtain enhance- 
ment, we divide by the background signal defined to be 
the experimental yield for the smallest available system, 
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FIG. 4: (color online) ^/ snn = 200 GeV yields of hyper- 
ons dA/dy and d~E~/dy, on left, and antihyperons dA/dy and 
dS /dy, on right, normalized with, and as function of, A, rel- 
ative to these yields in pp reactions. The lines are as in Fig.|2] 



The STAR collaboration has stated the cen- 



e.g.,A = 2. 

tral rapidity yield of dA/dy, dA/dy, dEr/dy and dE7/dy, 
at ^siN = 200 GeV in this way, see Fig.H l|. The ex- 
perimental reduced yields vary for large A, with A. This 
behavior cannot be due to a reinterpretation of the en- 
hancement as being due to the pp yield suppression, for 
it occurs at relatively large A. 

We overlay, over these data, the unsmoothed results we 
obtained with SHM. To normalize the SHM predictions, 
we use the central values of pp yields stated in Ref . 0] : 



= 0.066 ± 0.006, d(E~+ 3 )/dy = 0.0036 ± 
0.88±0.09andS + /S-= 0.90±0.09. Thus, 



d{A + A)/dy 
0.0012, A/A = 

we also extend the lines from A = 6.3, the smallest value 
considered in peripheral Au-Au interactions, to A = 2, 
where they assume by definition the value unity. Since 
the chemical equilibrium SHM fit is not valid below A = 
20 we do not continue this line. The edge, near to A = 
20, is signaling the change in behavior of the SHM fit 



associated with a possible phase change 0, ^| . 

The agreement of the chemical non-equilibrium result 
with all four data sets is remarkable. The reader should 
note that the relative yields and ^-dependence of the 
A, A, S, S are correctly described by SHM model param- 
eters derived from study of other particles. This is a 
remarkable accomplishment of the SHM. This remark 
applies to all three models considered, considering the 
magnitude of the normalization uncertainty, indicated 
at A = 2, in Fig.0] — here we note that we cannot 
with required precision extract the absolute yields from 
Fig-El 1-3] j an d enter these in Fig. [31 considering the dou- 
ble logarithmic representation. 

In this paper, we have explored the strangeness content 
of the RHIC fireball at its hadronization point found in 
Ref. |4(. We showed, as function of A, the behavior of the 
relative pre-hadron formation yields of light and strange 
quarks. Within the context of chemical non-equilibrium, 
the strangeness content increases rapidly with A. The 
reaction volume scales with A 2 / 3, , however, the residual 
variation is characteristic of phase properties and expan- 
sion dynamics. 

The increase in s/q with A, Fig.^ is imaged in the rise 
of the specific per hadron yields of strange antibaryons, 
Fig.0 with centrality. The variation of the strangeness 
yield with the reaction volume (impact parameter) is fur- 
thermore the source of the rapid rise of multi-strange 
(anti)baryon yield enhancement seen in Fig.^J 

The absolute strange (anti)baryon rapidity yields, 
Fig-El predicted for -^/snn = 200 GeV, distinguish 
the hadronization models: the yields of dEr/dy and 

dE + /dy are, for the chemical non-equilibrium case, 80% 
above those expected in the chemical equilibrium ap- 
proach. The high experimental hyperon yields reported 
at ^/snn = 130 GeV favor the chemical non-equilibrium 
hadronization. Future experimental results can as we 
have shown provide decisive distinction between the 
hadronization pictures considered. 
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